Kondo effect in a magnetic field and the magnetoresistivity of Kondo alloys 
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The effect of a magnetic field on the spectral density of a S = 1/2 Kondo impurity is investigated at 
zero and finite temperatures by using Wilson's numerical renormalization group method. A splitting 
of the total spectral density is found for fields larger than a critical value HdT = 0) ~ O.ST/f , where 
Tk is the Kondo scale. The splitting correlates with a peak in the magnetoresistivity of dilute 
magnetic alloys which we calculate and compare with the experiments on CexLai_xAl2, x = 0.0063. 
The linear magnetoconductance of quantum dots exhibiting the Kondo effect is also calculated. 
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The Kondo effect of a magnetic impurity interacting 
antiferromagnetically with conduction electrons has been 
the subject of theoretical and experimental investigations 
for several decades . It is most clearly manifested in the 
resistivity of electrons scattering from magnetic impuri- 
ties in dilute magnetic alloys, such as CexLai_xAl2, x = 
0.0063 ||], in which the resistivity shows an anomalous 
increase with decreasing temperature below a very low 
temperature Tk, the Kondo temperature, which is a 
scale determined by the many-body interactions. The 
phenomenon has recently acquired renewed interest in 
several contexts including transport through strongly in- 
teracting quantum dots and in the context of the 
dynamical mean field theory of strongly correlated lat- 
tice models Q. Although the thermodynamics of the 
Kondo model is well understood from exact solutions 
via the Bethe-Ansatz (BA) Q, spectral and transport 
properties have proven more difficult to calculate reli- 
ably, especially in the crossover region T « Tk where, 
neither perturbation theory about the Fermi liquid fixed 
point at T <C Tk, nor perturbation theory about the free 
fixed point at T 3> Tk, is adequate. A non-perturbative 
approach, such as the numerical renormalization group 
(NRG) method jo), is required to describe this crossover 
region. Recent developments in the NRG have yielded 
the zero field spectral and transport properties of the 
Anderson model in the Kondo regime 0,^. However, at 
finite magnetic fields, there are still very few calculations 
of these properties for a Kondo impurity . In this 
paper we calculate the field and temperature dependence 
of the Kondo resonance and the magnetoresistivity of a 
S = 1/2 Kondo impurity directly from the Kondo model 
by using the NRG method. The results are compared 
with available experimental data on dilute magnetic al- 
loys. We also calculate the equilibrium magnetoconduc- 
tance of a quantum dot in the Kondo regime. 

Model — Our starting point is the S = 1/2 Kondo 
model, which, for the purposes of this section, we write 
in the more general form of an anisotropic exchange fO] 
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The first term represents non-interacting conduction 
electrons and the second and third terms represent an 
exchange interaction between a localized spin 1/2 and 
the conduction electrons with strength J_l, J||. The mag- 
netic field, H, is taken to couple only to the impu- 
rity spin |]r^ . In the following we set g = /is = 
fcs = 1. The magnetoresistivity can be obtained from 
the T-matrix for electrons of spin a scattering from a 
Kondo impurity This is defined by the identity j^] 

Gk.k'A^) = 5kMQlk'A^) + Qlk,»%{^)Ql',k'»^ 

where Qk,k' ,cr{^) — {{ck,a',c^li o-)) is the retarded conduc- 
tion electron Green function and j^i the corre- 
sponding unperturbed Green function. The T-matrix is 
local due the local nature of the exchange interaction in 
(^. The impurity spectral density, A{lj, T, H), giving the 
Kondo resonance is the sum of up spin and down spin im- 
purity spectral densities, A{ui,T,H) = Ao-(ci;, T, iJ), 
where 

A^{Lj,T,H) = --Im%{uj + i5,T,H). (2) 

TT 

The linear response magnetoresistivity, p{T, H), due to a 
small concentration, c <C 1, of Kondo impurities is given 
by 

P-\T,H) . £yM.,T,H) , (3) 

where {uj,T, H) = cA„{llj,T, H) is the transport time 
of electrons of spin a, f is the Fermi function and m, n 
and e are the mass, concentration and charge, respec- 
tively, of the conduction electrons. An expression for the 
T-matrix is easily derived by considering the equation of 
motion of the Green function Gk,k',a- The result is 

^.H = Y(&) + ((o.;Oi)), 
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where co^o- = J^k^k.a-- This expression is also useful for 
calculating the properties of the Kondo lattice model via 
the dynamical mean field theory. 

Knowledge of the T-matrix also allows the magneto- 
conductance of a quantum dot in the Kondo regime to 
be calculated. The current, /, through such a quan- 
tum dot attached to two leads with chemical poten- 
tials ±6^/2 in the presence of a magnetic field H is 
given by I{T, H,V) ^ f- (w, T, H, V)if{u; + 

eV/2) - f{uj - eV/2)), where A^^^uj, T, H, V) is the spec- 
tral density of the dot generalized to a non-equilibrium 
situation V > JTst . A non-equilibrium approach is 
required to calculate A'^'^'^{uj,T, H,V) In experi- 

ments on quantum dots it is, however, possible to ac- 
cess the linear regime where jelV' <C Tk,H ^ and the 
relevant quantity is then the linear magnetoconductance 
G(T, H) = G{T, H,V = Q) = hmv^o dl/dV, 
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NRG calculations for G(T, H ^ 0) have been compared 
with experimental results on quantum dots in the Kondo 
regime and good agreement has been found Here, 
we present the field dependence of G{T, H) in the strong 
correlation limit Kondo regime of a quantum dot. 

Method — The technique we use to calculate 
Aa{uj,T,H) is Wilson's NRG method H extended 
to finite temperature dynamics 0. This is a non- 
perturbative method which gives the excitations and 
eigenstates of H at a decreasing set of energy scales 
uiN = A~^,{N = 1,2,..., A > 1), logarithmically 
spaced about the Fermi level = 0. It thereby al- 
lows a direct calculation of dynamic quantities, such 
as AaiujjT^ H), at a corresponding set of frequencies, 
uj ~ lon, and temperatures, ksT = Tn ~ i-^n, via the 
Lehmann representation. For further details see Ac- 
curacy is discussed below. 

Calculations — We now restrict ourselves to the case 
of relevance to Kondo alloys and quantum dots, J = 
J_L = ^ D, with D — 1 the half-bandwidth of the 
unperturbed conduction electron density of states per 
spin, N{uj), which we take to be a semi-elliptic. The 
number of states retained per NRG iteration was 462 
and a discretization parameter of A = 1.5 was used. 
We used J/D = 0.3, N{0)J = 0.19 for all calcula- 
tions, giving a Kondo scale of Tk — 0.0072, as deduced 
from the half-width at half-maximum (HWHM) of the 
Kondo resonance at T = (Fig. |l|). We shall use this 
definition of Tk throughout this paper. As expected, 
it is of the same order of magnitude as the expression 
= D{N{0)J)^/^e~^/^^°')-^ « 0.0023, which is valid for 
weak coupling A^(0) J <C 1. 

T-dependence of the Kondo resonance — The temper- 
ature dependence of the Kondo resonance in zero mag- 
netic field is shown in Fig. |l]. This is also known via NRG 




FIG. 1. Temperature dependence of the Kondo resonance 
a.t H — calculated from Eq. (P). The inset shows the peak 
height as a function of temperature. 



calculations on the Anderson impurity model |Q . By car- 
rying out calculations on the Kondo model we can focus 
more explicitly on the many-body aspects of this tem- 
perature dependence. The inset shows that the height 
of the resonance, which initially falls slowly at T ^ 
(as v4(0,0,0)(l — a{T /TkY)), decreases approximately 
logarithmically in a small range O.ITr- < T < Tr- close 
to Tk, with Tk defined from the HWHM of the T = 
resonance. The resonance falls to half its T = peak 
height at T « 0.5Tk. The u = T = value satisfies, for 
both zero and finite magnetic fields (see Fig. || below), 
the Fermi liquid relation (Friedel sum rule) 



A,iLo^O,T = 0,H) 
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where 5a{H) is the field dependent phase shift for elec- 
trons scattering from the impurity. This can be deduced 
from the Anderson impurity model ||l5[] , of which the 
Kondo model is the low energy part, or, by assuming a lo- 
cal Fermi liquid description for the Kondo model [|l6| . It 
provides a useful test on the accuracy of the calculations. 
On noting that the zero field phase shift in the Kondo 
model is 7r/2, we have A(0,0,0) = ^^^,,(0,0,0) = 
I/tt = 0.3183. The numerical result gives 0.3109 result- 
ing in a relative error of 2%. 

H-dependence of the Kondo resonance — The Kondo 
resonance "splits" in a large enough magnetic field — 
in the sense that the curvature in the total impurity 
spectral density at to — changes sign above a cer- 
tain field, see Fig. ||. Fig. || shows the evolution of 
the up spin part A^{uj,T, H) at T = as a function 
of field (A|(w,r,i7) = A^ {-uj,T,H)). The value at the 
Fermi level, given by Eq. (^), probes the field depen- 
dent phase shift 6-^ (H) . This has been calculated exactly 
from the BA solution of the Kondo model |^ by An- 
drei 0, via its relation, S^{H) = f (1 - 2M,(iJ)), to 
the impurity magnetization Mi{H). The inset to Fig. 
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FIG. 2. The up spin component of the Kondo resonance 
at T = for several magnetic fields. The inset compares the 
NRG and BA (BA) values for A^{lo = T ^ 0,H) for finite 
fields. 
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FIG. 3. Splitting of the Kondo resonance with increasing 
magnetic field at T = and T/Tk = 0.36. 

shows very good agreement between the NRG and BA 
for = 0, T = 0, H) for fields up to 5Tk. The same 

2% known error in the NRG result for the Friedel sum 
rule at H — was assumed for all fields and corrected 
for in the comparison. We estimate the shift, Eq{H), of 
the maximum in A-^iojjT = 0,H) to be Eo{H) = H at 
H -^Tk and EoiH) « 0.85F ± O.IH for H >Tk. 

The Kondo resonance depends strongly on tempera- 
ture and field, so it is not a priori clear that the total 
spectral density will split in any magnetic field either 
at T = or at higher temperatures. This depends on 
the competition between the reduction in height of the 
up spin and down spin components and their outward 
shift as a function of field. That this happens above a 
critical field. He, can be seen from Fig. ^ for two cases: 
T = and T/Tk = 0.36. Experimental evidence for 
this splitting comes from the temperature dependence 
of the magnetoresistivity of Kondo impurities, which we 
shall discuss below. Point contact spectroscopy on di- 



lute Kondo alloys ||T^ also indicates this splitting. At 
very low temperatures T < 0.25Tk we find that He is 
close to its T = value, 0.5Tk, and at higher tem- 
peratures is linear in temperature {Hc{T) « ST). We 
can compare these results with those obtained from a 
recent calculation of the one-spinon density of excita- 
tions for the Kondo model at T = 0. The low 
energy behaviour of this quantity should be close to 
the true spectral density Acr{uj 0,T = 0,H). For 
Hc(T = 0), which is a low energy property, the BA pc| ] 
gives giiBHc{T = 0) = 0.4858A:b7V in excellent agree- 
ment with our value of O-bknTK- Both approaches give a 
peak position Eq{H) lying close to but always below H . 

Magnetoresistivity — Fig. ^a, shows the temperature 
dependence of the magnetoresistivity together with a 
comparison to experimental results on the dilute Kondo 
alloy GexLai_xAl2, a; = 0.0063 We first describe the 
theoretical results. 

The magnetoresistivity for H = shows the behaviour 
known from calculations on the Anderson model 0|, in- 
cluding the T^ Fermi liquid corrections to the resistivity 
from its unitarity limit at T <C Tk Jl6t and the approx- 
imate logarithmic increase with decreasing temperature 
in a small range below T « Tk- On increasing the mag- 
netic field, a finite temperature peak appears at T « _ff 
in p(T, H) for fields larger than O.STr-. This reflects the 
splitting of the Kondo resonance, found above, for fields 
larger than 0.5Tr- and is physically reasonable, since the 
resistivity from Eq. (^) can be thought of as a broadened 
measure of the total spectral density. This statement can 
be made precise in the case of the conductance through a 
quantum dot in the Kondo regime, since from Eq.(^ we 
see that this quantity measures the field and temperature 
dependence of the Kondo resonance broadened by the 
derivative of the Fermi function. In Fig. ^ we quantify 
the difference between the magnetoresistivity of a Kondo 
impurity and the magnetoconductance through such an 
impurity by comparing them. From Eq.(|3|-^ we see that 
they are equal at T = O and all fields. They differ slightly 
only for fields comparable to or larger than Tk, and then 
only in the temperature range around T H. Hence, 
both quantities can be thought of as a "spectroscopic" 
measurement of the Kondo resonance |l3| . 

The dilute mag netic alloy CexLai^^Ala, x = 0.0063 |] 
was the best candidate we could find for a comparison of 
the magnetoresistivity to experimental results. It has a 
low lying doublet which can be modelled by a S" = 1/2 
Kondo model and both the normal part of the magne- 
toresistivity and the part due to disorder could be sub- 
tracted out 1^. The resulting magnetoresistivity at the 
lowest temperature was compared with T = BA |l^] 
results in JTot and found to be in good agreement, ex- 
cept at fields very small and very large relative to Tk- 
The finite temperature comparisons in Fig. ^ also show 
good agreement with the experimental results, except, as 
in the T = comparisons, at fields very small and very 
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FIG. 4. (a) Comparison of the calculated magnetoresistiv- 
ity curves (lines) and those of CexLai_xAl2, x — 0.0063 (sym- 
bols) Fields are in units of the theoretical Tk and corre- 
spond to H/kOe : 0, 1, 5, 10, 20, 50 of |] so our Tk = 0.8K 
(Tk ~ IK is cited in the experiment), (b) Comparison of 
the magnetoresistivity and magnetoconductance through a 
Kondo impurity (quantum dot). All curves in (a) and (b) 
were scaled by their unitarity limits (p(0, 0) and G(0, 0)). 



large relative to Tk (see below). There are no free pa- 
rameters in these comparisons. The curves were scaled 
by the unitarity value p(0, 0) and the relation between 
scales in experiment and theory was fixed by compar- 
ing the theoretical p{T = 0, H) versus H curve to the 
corresponding experimental one as in the comparisons 
to the T — BA calculations. This scaling factor for 
fields applies to all energies and was used to rescale the 
experimental temperatures onto the theoretical ones in 
units of our Tk- The large field dependence observed in 
the experiments between H — and H = IkOe cases 
(at T ^ Tk) could be due to the finite field extrapo- 
lation used to obtain the H — case (required because 
CexLai_xAl2, a; = 0.0063 is superconducting at _ff = 
[^). The experiments also give too strong a suppression 
of the magnetoresistivity with increasing field at large 
fields as compared with both our calculations and with 
the T — BA results The origin of this is unclear, 
but it may indicate that the real system is more compli- 
cated than an effective 5 = 1/2 Kondo model at these 
large fields. 

In summary, we have studied the splitting of the Kondo 
resonance in a magnetic field by calculating the T-matrix 
of the S = 1/2 Kondo model and we have shown that 
the splitting (in the total spectral density) occurs above 
a critical field Hc{T = 0) « Q.^Tk- This is manifested 
in a peak at finite temperature for H > 0.5Tk in the 
magnetoresistivity of Kondo alloys and in the magne- 
toconductance of a quantum dot in the Kondo regime. 
Both quantities give "spectroscopic" information on the 
Kondo resonance and we showed explicitly that they are 
almost equal for most fields and temperatures. Quanti- 



tatively good agreement with the magnetoresistivity of 
the Kondo alloy CexLai_xAl2 was found with no free pa- 
rameters. Our results should be of use in interpreting 
experiments on the field dependence of the conductance 
of quantum dots in the Kondo regime. 
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